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In experiments conducted recently at MIT on Na Bose- 
Einstein condensates [S. Inouye et al, Nature 392, 151 (1998); 
J. Stenger et al, Phys. Rev. Lett. 82, 2422 (1999)], large loss 
rates were observed when a time-varying magnetic field was 
used to tune a molecular Feshbach resonance state near the 
state of a pair of atoms in the condensate. A collisional deac- 
tivation mechanism affecting a temporarily formed molecular 
condensate [see V. A. Yurovsky, A. Ben-Reuven, P. S. Juli- 
enne and C. J. Williams, Phys. Rev. A 60, R765 (1999)], 
studied here in more detail, accounts for the results of the 
slow-sweep experiments. A best fit to the MIT data yields 
a rate coefficient for deactivating atom-molecule collisions of 
1.6 x 10~ 10 cm 3 /s. In the case of the fast sweep experiment, 
a study is carried out of the combined effect of two com- 
peting mechanisms, the three-atom (atom-molecule) or four- 
atom (molecule-molecule) collisional deactivation vs. a pro- 
cess of two-atom trap-state excitation by curve crossing [F. 
H. Mies, P. S. Julienne, and E. Tiesinga, Phys. Rev. A 61, 
022721 (2000)]. It is shown that both mechanisms contribute 
to the loss comparably and nonadditively. 
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I. INTRODUCTION 

Most properties of a Bose-Einstein condensate (BEC) 
are determined by interatomic interactions (see Refs. 
jl],D). These interactions are responsible for the char- 
acteristic nonlinear term in the equations of motion of 
the condensate, whose magnitude depends on the colli- 
sional elastic scattering length. Recent experiments [||^] 
on optically-trapped BECs drew attention to the effects 
of Feshbach resonances on the properties of the conden- 
sate, as scattering lengths are strongly modified by the 
presence of a resonance. More particularly, these experi- 
ments sought to modify these effects by application of a 
magnetic field (see Refs. One of the rather aston- 

ishing results observed was a dramatic loss of condensate 
population as the magnetic field was varied so that the 
ensuing Zeeman shift made the system pass through a 
resonance, or approach it closely. 

A Feshbach resonance may exist when the energy of 
a pair of atoms in the condensate is close to that of a 
metastable molecular state Na2 (to). Then the scatter- 
ing length varies strongly as a function of the energy 



mismatch between the two states. This mismatch can 
be controlled by applying a varying magnetic field. The 
energies of the two states can be brought closer to each 
other, as the two states have different Zeeman shifts. The 
MIT experiment (HJJ] strived to study the effect of the 
Zeeman shift by applying a time-varying magnetic field 
in two distinct procedures: (a) a fast sweep through the 
resonance, using a fast ramp speed of the magnetic field, 
and (b) a slow sweep, using much lower speeds, in which 
the ramp is stopped short of crossing the resonance. The 
latter procedure is then repeated by using different val- 
ues of the stopping value of the magnetic field, and is 
carried out on both sides of the resonance. Both types 
of experiment resulted in a large condensate population 
loss. 

This work is devoted to the study of possible mecha- 
nisms leading to this loss. Preliminary results were pre- 
sented in Refs. [0J|], suggesting the mechanism of colli- 
sional deactivation. Another mechanism, involving trap 
excitation by a two-step curve-crossing process, was sug- 
gested in Refs. [6 10 1 . We study below in more detail the 
combined effect of both mechanisms. As the MIT experi- 
ments |^ were conducted with Na atoms, we shall refer 
here for definiteness to Na only, though this study may 
be extended to similar systems. 

Given a stationary Zeeman shift, a population of 
molecules can be formed as a temporary stage in the 
elastic process 



Na(BEC) 



f- Na(BEC) 
Na(BEC) 4 



-> Na 2 (to) 
Na(BEC). 



(1) 



However, in the absence of other intervening interactions, 
or of a time- varying field, this molecular population can- 
not persist, and no loss would occur. 

The first loss mechanism considered here involves the 
deactivation of the resonance state, which is usually a 
highly excited vibrational level in a given spin state of 
the pair, by an exoergic collision with a third atom of 
the condensate [§|,|],D, 



Na 2 (m) + Na(BEC) -> Na 2 (d) + Na(hot), 



(2) 



bringing the molecule down to a stable state d, and re- 
leasing kinetic energy to the relative motion of the re- 
action products. Although the collision occurs with a 
vanishingly small kinetic energy, rates of such inelastic 
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processes remain finite at near-zero energies p). This 
process naturally depends on the initial density in the 
condensate. A variant of this process, involving deacti- 
vation by collisions with another molecule (rather than 
an atom) , of the type 



Na 2 (to) + Na 2 (m) -> Na 2 (d) + Na 2 (it) 



(3) 



would require a significant molecular density. The two 
molecules emerge in two states d and u, where u can be 
a stable molecular state above d, or a continuum state 
of a dissociating molecule. The reaction can take place 
as long as the corresponding internal energies obey the 
inequality + E u < 0, where the internal energy of to 
serves as the zero reference point on the energy scale. 
A particularly effective reaction of type (^) would occur 
in the near-resonant case, in which < E u < \Ed\- A 
typical example, common in VV-relaxation, is that of v + 
v — ► (v + l) + (v — 1) (where v is the vibrational quantum 
number of the state m). In this example, the kinetic 
energy is provided by the vibrational anharmonicity. 

Both reactions (^) and (|3j) are thus exoergic, provid- 
ing products with sufficient kinetic energy to escape the 
trap, as the characteristic transition energies exceed the 
trap depth. The kinetic energy may even be sufficient to 
produce an additional loss mechanism — secondary col- 
lisions of the reaction products with condensate atoms 
(see Ref. ||). The loss mechanism described here can be 
enhanced by bringing close to each other the energies of 
the two states involved in reaction (|l|) — the BEC state 
of an atom pair and the resonant molecular state to — 
by the application of a time- varying Zeeman shift. It 
is not necessary that the energies of these states should 
cross. An actual crossing of the two states can cause an 
irreversible transfer of population from the condensate to 
the molecular states (see Refs. ^-|l0|). This crossing is 
also necessary, as the first step, in the other mechanism 
referred to earlier — that of excitation of the trap states 
by a two-step crossing. The first crossing between the 
condensate and molecular states can occur in two direc- 
tions, either by letting the molecular state move upwards 
in its energy, with respect to the condensate atom-pair 
state, or by letting it move downwards. In both cases 
the loss of the molecules can proceed via the deactiva- 
tion mechanism. But only the upwards move alone can 
initiate the excitation mechanism. The second crossing, 
at a higher energy, then causes transfer of population to 
higher atomic states — bound and continuum trap states 
in a condensate imbedded in an optical trap ]Tc| ] , or con- 
tinuum states in a free condensate ||. This process is 
accompanied by an increase in energy 



Na(BEC) + Na(BEC) Na 2 (to) -> Na(hot) 



Na(hot). 

(4) 



the first one along the time axis in a Z-shaped forma- 
tion jl2|]. However, this effect is negligibly small in the 
present case.) 

It should be made clear that both crossing directions 
have been taken into account in Ref. |lQ] . However, that 
work does not specify what happens to the molecular 
population in the downward move. In principle, once the 
Zeeman shift undergoes the first crossing, the formation 
of a molecular population in the resonant state can be 
considered as a valid loss channel, no matter whether 
there are deactivating collisions or not. In that case, the 
loss rate should be independent of the deactivation rate. 
Our analysis here shows that this is generally not the 
case. In the case of the fast-sweep experiment, the loss 
would become independent of the deactivati on rat e only 
under certain conditions (referred to in Sec. II C 2 below 
as the "asymptotic" conditions). 

It is now understood |^|| that in the slow-sweep ex- 
periment the loss is produced almost exclusively by deac- 
tivating collisions, such as the atom-molecular collisions 
described by (||). We study here also the added effect 
of molecule-molecule collisions described by (|^) . In the 
case of the fast-sweep experiment, it has been claimed 
P,|To[ that the main cause of loss is due to the excitation 
processes of the kind described by (Q) , but it was shown 
H that deactivating collisions cannot be discounted as a 
contributing mechanism. 

This paper therefore aims to study the effect of combin- 
ing the two kinds of mechanisms (condensate excitation 
vs. deactivating collisions) together. One of the major 
conclusions (discussed in Sec. IV below) is that the two 
processes may actually compete nonadditively with each 
other, rather then contribute additively to the loss pro- 
cess. 

The paper begins, in Sec. [n| by an expansion of the 
theoretical analysis used in ||] to describe the effect of de- 
activating collisions. We show, among other things, how 
the coupled equations of the Gross-Pitaevskii type for the 
atomic and molecular condensates, introduced and stud- 
ied earlier by Timmermans et al. Q, can be derived by 
an elimination of the product states of the deactivation 
process (the so-called "dump" states) . The equations are 
then extended to include molecule-molecule collisions of 
type (|). In Sec. fil we add to the deactivation model an 



(To be more precise, two-step excitation can in principle 
occur also in a downward move by the so-called "counter- 
intuitive" process in which the second crossing precedes 



effect representing the outcome of the excitation process. 
The results are presented and discussed in Sec. IV, in 
comparison with the MIT experiments. 



II. THE DEACTIVATION MODEL 

A. Hamiltonian and variational procedure 

Let us consider an optically-trapped BEC exposed to 
an external homogeneous time-dependent magnetic field 
B (t) used to tune a vibrationally excited molecular state 
to to a Feshbach resonance with the state of a pair of 
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unbound condensate atoms. In order to write down an 
hamiltonian for such a system, including the molecular 
dump states, we must introduce field annihilation op- 
erators of the atoms V>(r), of the molecular resonant 
state tp^ (r m ), and of the lower and upper dump states, 
( rm ) an d K ( r m), respectively [see discussion follow- 
ing Eq. ([}])]. The hamiltonian can then be written as 



H 



d 3 r^ + (r) H a <>p (r) + V el 

/ d 3 r m ^2 i>t ( r "0 &ai>a (r m ) 

« ™ A 



a. — m , u . d 



+v h + v+ + 



E(^ + ^ + )+E(^ + ^)- ( 5 ) 



The terms 



H a = t^P 2 + V a (r) - p a B (t) 
Am 



H a = —pi + V a (r m ) - p a B (t) , 
4m 



(6) 



(where a = m, u, or d, includes the resonant state) 
are the hamiltonians for the noninteracting atoms and 
molecules. Here V a (r) and V a (r m ) are the correspond- 
ing atomic and molecular energies as functions of the 
position r of the atomic (or r m of the molecular) center 
of mass, whose values include the optical trap potentials 
and the position-independent differences of internal en- 
ergies in the absence of the trap. Also, /i a and p a are the 
corresponding magnetic moments. 

Since the atoms and molecules are treated here as in- 
dependent particles, the interaction responsible for the 
atom-molecule coupling [reaction ([!])] can be written in 
the general form 



Vk 



d 3 rd 3 r'V h (r - r') 



^(r)^(r') 



(7a) 



in which the molecule, created as an independent parti- 
cle preserves the position of the center of mass. How- 
ever, considering that the interaction is localized within 
a range of atomic size, negligibly small compared to the 
condensate size and the relevant de Broglie wavelengths, 
one can use the approximation of zero-range interaction 
Vh (r — r') = gS (r — r'), and represent the interaction in 
the simpler form 



V h = g d 3 r^ (r) ^ (r) i> (r) . 



(7b) 



The same arguments are applicable to the terms in Eq. 
(||) representing the deactivating collisions Vd and V u d 
[reactions (||) and (|3|), respectively]. However, the use of 
a zero-range interaction would lead to a divergence in the 



ensuing calculations [see discussion following Eq. (g(j) be- 
low] . Therefore we keep these interactions as finite-range 
functions of the distance between the reaction products, 
writing 



V d = I d 3 rd 3 r m d d (|r - r m |) V + (r) (r m ) 
r + 2r m \ - / r + 2r„ 



V ud 



3 J rm V 3 
d 3 ri d 3 r 2 d u d (|n - r 2 |) (n)^+ (r 2 ) 
n + r 2 \ 7 fr 1 +r 2 



(8) 



(9) 



Here, as in (7a) the position of the center of mass is pre- 
served, but the finite-range nature of the interactions is 
retained in the functions dd (p) and d u d (p), whose actual 
shape will be discussed further down. 

Finally, the part of the hamiltonian associated with 
elastic collisions (see Ref. f|), 



V, 



d 3 r 



-y ^ + (r) V> + (r) ip (r) tp (r) 
. U r . 



E ^« w ( r ) ( r ) ^ a ( r ) 



+U am E V> (r) i>+ (r) i, a (r) V> (r) 



. (io) 



includes terms proportional to the zero-momentum atom- 
atom, molecule-molecule, and atom-molecule interac- 
tions, 



U a 



U„ 



U a 



(ii) 



where a a , a m , and a am are the corresponding elastic scat- 
tering lengths, and the different numerical factors in the 
numerators reflect the different reduced masses. 

We outline here the derivation of mean- field equations, 
involving c-number fields, that represent all the actively 
participating states listed above. This is accomplished 
by an extension of a well-known variational method (see 
Refs. |i~3 14|). Let us introduce the trial wavefunction 



1+ / d 3 rd 3 r m ^ip d {Y,Y m ,t)'4> + (r)V>+ (r m ) 

d 

+ I d 3 r 1 d 3 r 2 ^^ ud {Y 1 ,Y 2l t)i ! + (n)^+ (r 2 ) 



x|</?o, <p m )- 
The factor 



tid 



(12) 



Wo, 9m) 



exp 



d 3 r[w (r,t)^+ (r) 

+p m (r,t)^+(r)]}|0) (13) 
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is a coherent state, formed by a product of exponential operators, involving atomic ((fin) and molecular (ip m ) condensate 
states, and operating on the vacuum state |0). The linear factor preceding it in Eq. ( |l2|) includes the fields ifid (r, r 2 , t) 
and (fi u d (ri, r 2 , i) which are the correlated states of the products in reactions (||) and (||), respectively. 

The linear form of the latter factor forces the trial wavefunction ((l^) to take into account only one-particle occupation 
of the non- resonant molecular states u and d, as a constraint. This approximation is based on the assumption of fast 
removal of "hot" particles from the trap. In contrast, the occupation of the resonant state (to) is allowed to reach the 
order of magnitude of the condensate-state occupation (as our calculations verify). Therefore if m describes a coherent 
molecular condensate. 

Another constraint, that follows from the large energy difference between the dump states and the resonant state, 
is the condition 

[ d 3 r<fiZ(r,t)ip d (r,r m ,t) =0 (14) 

[see discussion followi ng Eq . ( |2(i| ) for justification] . The trial wavefunction ( |l2| ) can now be substituted into a variational 
functional (see Refs. |l^ , |l4| ) 

oo 

d V 1 * i / , 15 
($1$) 

Neglecting terms of the order of tp^ and tp^, and taking (|lj) into account, the use of a standard variational procedure 
(see Ref. jl4j) then leads to a set of coupled equations for the atomic (<po) an d molecular (<p m ) condensate fields (or 
"wavefunctions"), as well as for the dump states (ifid and <pud)'- 

ih(fi (r, t) = (ll a + U a \(fi (r, t) | 2 + U am \ifi m (r, t) | 2 ) <p (r, t) + 2g*tp* (r, f) ip m (r, t) + Q (r, t) tp* m (r, t) (16a) 

i?iy> m (r,i) = (H m + U am \tp (r,i) | 2 + U m \tp m (r,t) | 2 ) tp m (r,i) + gp% (r,t) 

+Q (r, t) p* (r, i) + Q m (r, t) ^ (r, t) (16b) 

— " 2 + — 
2rrJ >1 4to 

— " 2 + — 
Am^ 1 4to 



ify&d (ri,r 2 ,t) = I -^-p 2 + ^—p\ - E d ) ifid (ri,r 2 ,t) + d d (|n - r 2 |)^o j ri H ~ 2r2 ,t ) p m | r i + 2r 2 ^ i 



itufi ud (ri,r 2 ,t) = [ -f-ffi + ^—pl - E u — E d ) (fi u d (ri,r 2 ,t) + d ud (|n -r 2 |)p£, I ^— ) , (16d) 



where 

ri + 2r 2 



Q(r,t) = J d 3 r 1 d 3 r 2 ^d* d (\v 1 ~ v 2 \) i Pd (r 1 ,r 2 ,t) 8 (r 



ri + r 2 



(17) 
(18) 



Q m (r,i)=2 / d 3 r 1 d 3 r 2 ^2d* ud (\ri~r 2 \)ifi* m (r,t)ip u d(r 1 ,r 2 ,t)5 (r - 

J ud ^ 

The term s corresponding to elastic collisions, dependent on \<fid\ 2 and |(/5 U( ;| 2 , which should appear in E qs. (16a) 
and (16b), are omitted since they are negligible compared to the terms including ipn and ip m . In Eqs. (16c) and 
( |l6cj ) the terms corresponding to elastic collisions, Zeeman shifts, and trap potentials are neglected compared to the 
position- independent internal energies, denoted here Ed and E u . 



B. Dump state elimination 



The procedure used to eliminate the dump s tate s is similar to the Weisskopf-Wigner method of the theory of 
spontaneous emission (see Ref. [Q). Equation (16c) is of the form of a Schrodinger equation for two free particles 
with a source (the last term in the right-hand side). Such an equation can be solved by applying the Green's function 
method for free particles, with the result 



ifid (ri,r 2 ,£) = - 



(2tt) d n 



x exp 



zK 



dt' / d 3 Kd 3 kexp 



+ ik (r x - r 2 ) 



ffiK 2 



ri + 2r 2 



V 6to 



ihk 2 

Am 



Ed 



(t-f) 



d 3 pd d (p) e~ ik " / d 3 Re- iKR <fi m (R, t') p (R, f ) . 



(19) 
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Here R is the center-of-mass position of the three-atom system, p is the radius vector of the reaction products, and 
K, k are the corresponding wavevectors. Since ip m (R, i) and tpa (R, t) are condensate wavefunctions, the Fourier 
transform of their product [the integral over R in Eq. (|l9|)] vanishes if K > 1/6, where b is a characteristic size of the 
condensate. Therefore Ti 2 K 2 / (6m) < haj tTa , p is negligible compared to Ed 3> ftwtrap, where Wtrap is the trap frequency. 



This fact allows us also to neglect the time dependence of ip 7 
obtain the simplified expression 



(R, t) and ipo (R, t) in the integration over t', and thus 



(fid (ri,r 2 ,f) 



1 /r! + 2r 2 

rVo 5 ,* if r , 



ri + 2r 2 



(2tt) 3> ' u V 3 

The atom-molecule pair is thus formed with a momen- 
tum hkd — J AmEd/2) of relative motion. The function 
ipd (ri , r 2 , t) is a rapidly oscillating function of the coor- 
dinates, and therefore condition ( ]14|) is justified. 

Substituting Eq. @ into Eqs. l|7|), and @ and in- 
troducing a Fourier transform of the function dd {p) 



dd (k) = I d 3 pd d (p) exp (-ikp) , 

one obtains 

Q (r, t) = - (S + ihj) ip (r, t) ip m (r, t) , 



where 



|d„ (k) 



3h 2 k 2 /(4m)-E d -iO' 



(21) 



(22) 



(23) 



^kdldd (k d ) 



2m 



OO 



k 2 \d d {k)\ 
k 2 -k 2 



A similar analysis, starting from Eq. (16d), gives 
Q m (r, t) = - (S m + ifij m ) ip 2 m (r, t) , 

where 



In 



— Vfc 



u . d 



2m 
^2 



{k U ci) 

fc 2 M„ d (fc)| 2 



(24) 
(25) 

(26) 

(27) 
(28) 



and k ud = \/2m ( E u + E d )/h Substituting Eqs. flgj ) 

and into Eqs. (16a) and ( 16b ) one finally obtains a 
pair of coupled Gross-Pitaevskii equations (see Ref. |16|) 



d 3 kd 3 p 



d d (p) exp [zk (r 2 - ri - p)] 
3?i 2 fc 2 /(4m) -Sd-iO 



(20) 



Using the well-known identity [x — iO) 1 = Vx~ x + 
iird (x) , where V denotes the Cauchy principal part of 
the integral, allows us to obtain explicit expressions for 
7 and 5: 



ihip = [fl a + U a \<Po\ 2 + U am \ip m \ 2 ) ipo + 2g*ip* Q (p m 

-(5 + ^ 7 )|^ m | 2 v3o (29a) 

|2 , tt I,. \2\ ,„ i „,„2 



ihlp m = [H m + U am \fo\ + U m \fm\ j <Pm + 9<Po 

- [(6 + ih-y) \<po\ 2 + {S m + ihjm) \tp m \ 2 ] <p m - (29b) 

The parameters (5, 7, S m , and 7™, which are expressed in 
terms of d d and d ud [see Eqs. @, ©, @, and (§§)], 
describe the shift and the width of the resonance due 
to the deactivating collisions with atoms and molecules, 
respectively. The parameters 7 and "f m are one half of the 
corresponding rate constants. Since the strengths of the 
deactivating interactions are unknown, the parameter 7 
will be extracted below from the experimental data, and 
7 m will be used below as an adjustable parameter. The 
shifts S and 6 m can be incorporated in the interactions 
Uam and U m , respectively. 

Equations ( p9| ) are similar to those presented recentl 
by Timmermans et al. Q. Among other things, Ref. 
shows that in the case of a time-independent magnetic 
field and large resonant detuning, neglecting the decay 
described by the imaginary terms, Eqs. (29) can be re- 



duced to a single Gross-Pitaevskii equation with an ef- 
fective scattering length a a [1 — A/ (B — B Q )], where the 
parameter A is related to the atom-molecule coupling 
constant g of Eq. (ffq ) as 



1.91 



27r?i 2 |a a |/LtA/m. 



(30) 



Values of A for Na were calculated in Refs. 
extracted from the experimental data in Refs. 



C. Density equations and approximate solutions 



Let us introduce the new real variables 



n (r, t) = \(p (r, t) | 2 , n m (r, t) = \ip m (r, t) f 
u(r,t) = 2Bfi (g^ 2 (r,t) </4 (r,t)) /h, 
v (r, t) = -21m (gip 2 (r, t) ip* m (r, t)) j%. 



(31) 



The time evolution of these variables can be described 
by a set of real equations, similar to the optical Bloch 
equations where n and n m act like "populations" and u 
and v like "coherences" . When the kinetic energy terms 
are neglected, in accordance with the Thomas-Fermi ap- 
proximation (see Refs. 0fl[), one obtains from Eqs. (|2T 
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n = 2v- 2T a n (32a) 

n m = —v - 2T rn n m (32b) 

v = Du- (2T a + T m ) v + 2\g\ 2 n (4n m - n) /h 2 (32c) 

u = -Dv - (2T a + T m ) u. (32d) 

Here 

D(r,t) = {V(r)-fiB (t) 

+2 [U a n (r, t) + (U am - 5) n m (r, <)] 

- [{U am ~S)n (r, t) + {U m - 6 m ) n m (r, t)]}/h, (33) 

V (r) = 2V a (r) -V m {r), fi = 2fi a - fi m , 

and 



T a (r, t) = jn m (r, t) , 

r m (r, t) = jn (r, t) + 7 m n m (r, t) . 



(34) 



In the Thomas-Fermi approximation the functions 
n(r,t), n m (r,t), v(r,t), and u (r, t) depend on r only 
parametrically. The set of four real equations ( p{2] ) 
can then be solved numerically using as initial condi- 
tions either an r-dependent (for example, a steady-state 
Thomas-Fermi) distribution, or an r-independent (homo- 
geneous) distribution equal to the mean trap density. 

Nevertheless, certain properties of the solutions can be 
derived analytically from Eqs. (|3^), without recourse to 
numerical solutions, whenever the following "fast decay" 
conditions hold: 

liB « Jir m (D + T 2 JD) , r m » r Q (35a) 
D 2 +T 2 m > 6\g\ 2 n/h 2 . (35b) 

These conditions mean that the relaxation of n m , v, and 
u is much faster compared to that of n and to the rate of 
change of the energy, caused by the magnetic field with a 
sweep rate B. Therefore the values of the fast variables 
can be related to a given n value, using a quasi-stationary 
approximation, by 



D 



-v, 



2\g\ 2 n 2 T m 
h 2 (D^+Tiy 



(36) 



\g\ 2 n 2 



h 2 (D 2 + VI 



and the condition ( 35b ) leads to n m <C n. As a result, 
a single non-linear rate equation for the atomic density 
can be extracted. When terms proportional to the atomic 
and molecular densities in D [see Eq. d^)] are neglected, 
the resulting rate equation is 



h(r,t) 



6|<?|V 3 (r,*) 



[V(r)-»B(t)Y + {n ln {v,t)y 



(37) 



(The neglected terms in D effectively add an extra shift 
to the resonance, but its contribution is hardly noticed 
in the present problem.) 



Equation (37) has a form analogous to the Breit- 
Wigner expression for resonant scattering in the limit of 
zero-momentum collisions (see Ref. |l(J). In the Breit- 
Wigner sense one can interpret 2hrjn as the width of the 
decay channel, while the width of the input channel is 
proportional to |g| 2 . This observation establishes a link 
between the macroscopic approach used here and micro- 
scopic approaches that treat the loss rate as a collision 
process. However, Eq. ( |37| ) differs from the usual Breit- 
Wigner expression by a factor | , associated with the loss 
of a third condensate atom in the reaction (0). 



1. Approaching the resonance 

Very close to resonance (e.g., for conditions prevail- 
ing in the experiment Q, where B (t) is within 1/iT of 
resonance) the behavior of Eq. (]37]) effectively attains a 
1-body form linear in n. But as long as we stay out of this 
narrow region, by obeying the " off-resonanc^ condition 
(which holds in the slow-sweep experiment [0fl]), 



?i 7 n(r,i)« \V{r)-»B (t) |, 



(38) 



we can write Eq. ( p7| ) (to a very good approximation) in 
the 3-body form n = —K3 (r, t) n 3 , where 



Ko = 



I2ir?i 2 \a a \-fk 



m[i [B (t) - V (r) /n] 



(39) 



The dependence of Eq. ( p9| ) on the scattering length a a 
follows from Eq. (^0|). A similar expression has been ob- 
tained in Ref. ||, but the loss of the third condensate 
atom in the deactivating reactio n (H ), described by the 



term proportional to r a in Eq. (32a), was neglected. In 
the fast-decay approximation, in which h m is diminish- 
ingly small, this neglected term is equal to —v. Added 
to the — 2v, this makes K% of Eq. (|39"|) larger by a fac- 



tor I than the corresponding expression in Ref. |6j . This 
omission has been corrected in Refs. J|||. 

When the magnetic field ramp is assumed to vary lin- 
early in time, starting at to and ending at t, and Eq. 
( p8| ) applies throughout the ramp motion (i.e., by avoid- 
ing passage through the resonance), the rate equation 
can be solved analytically. Using Eq. (B9h one then gets 



n(r,t) = n(r,t ) 



1 + 24nh 2 \a a \Ajn 2 (r,t ) 
1/2 



(t - t ) I (m[iB 2 tto 



(40) 



where B is the magnetic-field ramp speed and the extrap- 
olated time of reaching exact resonance is chosen to be 
t = 0, so that both t and to have the same sign. We shall 
refer to the combination of Eqs. (35) and (pq), that leads 
to conditions (|39"|), as the "three-body" approximation. 



G 



2. Passing through the resonance 

The three-body approximation of Eqs. (^) and (|39| ) 
does not hold very close to resonance, and is therefore 
inapplicable to the description of the fast-sweep experi- 
ment, in which the Zeeman shift is swept rapidly through 
the resonance, causing dramatic losses (see Refs. 
Nevertheless, the fast decay approximation may still 
be valid. A simple analytical expression can then be de- 
rived for the condensate loss on passage though the res- 
onance if, in addition, the magnetic field variation lasts 
long enough to reach the 11 asymptotic? condition 



(41) 



where SB is the total change in B accumulated over the 
sweep. This condition allows the extension of the ramp 
starting and stopping times to =Foo. 

The asymptotic behavior of n (t) and n (t) in the com- 
plex t plane, as \t\ — ► oo, is constrained by consistency 
requirements. Consider the four possible asymptotic re- 
lations between n and t shown in the first column of Table 
H Equation |37]) forces n to attain the form in the sec- 
ond column, from which it follows that n should attain 
the form in the third column. Obviously, cases (c) (for 
Ret > 0) and (d) (at all t values) are not self-consistent. 
Therefore, the asymptotic solution may attain only one 
of the forms complying with cases (a), (b), or (c) (the 
latter for Ret < only). 

One can now evaluate the variation of n (r, t) in the 
infinite time interval (— oo,oo). Let us rewrite Eq. ( pT^ ) 
in the form 



dn 
n 2 



6jn\g\ 



liBt) +(h-yn) : 



-dt, 



(42) 



where V (r) is removed by our choice of the origin on 
the time scale. We then integrate the left-hand side with 
respect to n from n (r, — oo) to n (r, oo) and the right- 
hand side with respect to t from — oo to oo, considering 
n as a well-defined function of t. The latter integral may 
be evaluated by using the residue theorem, closing the 
integration contour by an arc of infinite radius in the 
upper half-plane. The integral along this arc vanishes 
according to the asymptotic behavior of n considered in 
all self-consistent cases of Table [| 

The final result does not depend on 7, and on the self 
-consistent case studied, and has the form (valid for all 
positions r) 



n (r, 00) 



n (r, —00) 



1 + sn (r, —00) ' 

67r| ff | 2 _ 12ir 2 h\a a \ A 
hn\B\ m 



(43) 



\B\ 



The product sn in Eq. ( p3[ ) would be proportional to the 
Landau-Zcncr exponent for the transition between the 



condensate and the resonant molecular states whose en- 
ergies cross due to the time variation of the magnetic 
field, if one could keep the coupling strength gip con- 
stant. However, for the non-linear curve crossing prob- 
lem represented by Eq. (29), the Landau-Zener formula 
is replaced by Eq. (|43|), which predicts a lower crossing 
probability, since the coupling strength gipo decreases, 
along with the decrease of the condensate density during 
the process. 

The asymptotic result ( [43] ) describes the decay of the 
condensate density. Assuming a homogeneous initial 
density within the trap, Eq. (|4^) applies also to the loss 
of the total population N (t) = J n (r, t) d 3 r. 

An asymptotic expression for the total population can 
also be found when the homogeneous distribution is re- 
placed by the Thomas- Fermi one (see Ref. fL7|). In this 
case, given no is the maximum initial density in the cen- 
ter of the trap, one obtains 



1 



at (00) 15 r 

JV(-oo) ~ 2sn { 3 



1 

sn 



1 



2S71Q 



1 

sn a 



x In 



1 

sn 



+ 1 / 



1 

sn 



(44) 



III. INCLUSION OF LOSS BY EXCITATION 

The other mechanism of molecular condensate loss, 
considered in Ref. ]l(J, involves the decay of the res- 
onant molecular state by transferring atoms to excited 
(discrete) trap states or to higher-lying non-trapped (con- 
tinuum) states. This decay is possible when the potential 
energy of the resonant state V m — ^, m B (t) exceeds the en- 
ergy of the atom pair formed. A simpler version of this 
mechanism has been studied in Ref. ||. 

Let us consider, for a moment, the decay of the 
molecule as a "half collision" , ignoring the finite size of 
the trap. The excess energy E of the released pair, where 
E = E (t) = /.iB (t) — V at the time of release, determines 
a "half-width" of the Feshbach resonance. According to 
Eq. (27) of Ref. jio) this energy-dependent half-width is 
given by 



r cr (E) 



h 2 



(45) 



It should be recalled that here always E > 0. The cou- 
pling of the resonant molecular state with an excited 
state v of an atom pair in the trap can be described 
by a coefficient g v . This coefficient can be related to r cr 
through (see Ref. pi) 



2 h r 

9v = -r c 

7T 



de v 
dv ' 



(46) 



where e v is the pair excited-state energy measured from 
the ground trap state and de v /dv measures the distance 
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between the trap states in the vicinity of v (i. e., the 
inverse density of states) . 

As the magnetic field B (t) = Bq + Bt rises above the 
resonant value -Bo, a crossing starts to occur between 
the resonant molecular state and the excited trap levels. 
Neglecting molecular collisions and motion, the resonant- 
state wavefunction (p m (t) is propagated from ip m (to) ac- 
cording to the rule (see Ref. fig ]) 



Vm (t) = <p m (to) exp 



(47) 



Although g v , for a given v state, is time- independent, 
the sum taken over all v states, bounded by the interval 
[iBto < e v < \iBt, is time-dependent. Whenever the 
amplitude of each crossing is small, i.e., when 

\g v \ 2 / (Vb) « 1, (48) 
the sum in Eq. (f47|) can be replaced by an integral, giving 



^ m (t) = <p m (to) exp ( ~~^§J ^ de%v dt I ' (49) 

where (de v /dv) / fiB measures the time interval between 
sequential crossings. Differentiation of Eq. p9| ) with re- 
spect to t gives the following expression 



dlfirr, 

at 



\9v\ 



dv 
de v 



<Pr, 



-r cr (e«) fr, 



(50) 



It should be kept in mind that this result is valid only 
in cases in which the energy gap between the atomic and 
molecular states increases rather fast; i.e., when fiB > 
and the condition (fj§|) is obeyed. In the case of \iB < 0, 
transitions from the ground trap state to excited ones 
are counterintuitive (see Ref. fL2fl ), and become negligible 
when 



SB^> 



(51) 



trap 



where SB is the range of variation of B extended over 
both sides of the resonance. 

Thus, whenever the energy gap between the resonant 
molecular and condensate states is positive, it increases 
rather fast [following Eq. (4§)], and Eq. ( |5l"| ) is obeyed, 
one can account for the decay of the resonant molec- 
ular state into excited trap states by adding a term 
—ihT c 



(jxBtj ip m to the right hand side of Eq. ( 29b ) or 
by substituting T m = jn + "f m n m + T CI (idBt) in Eqs. 



(p2|). In this way, the two mechanisms can be combined 
in a single formalism. Calculations made using this for- 
malism are discussed in the following section. 



IV. RESULTS AND DISCUSSION 

Calculations have been carried out on the loss of atoms 
for the case of the Na BEC experiments p,H] using both 




and (jWj)] and 

The parameters used 



analytical results [Eqs. (|3£ 
numerical solutions of Eq. flg 

to describe the system, reported previously 10|9|], are 
a a = 3.4nm, /i = 2/i a — (j, m = 3.65/xb (where fiB — 
9.27 x 1(T 24 J/T is the Bohr magneton), and A = 0.95//T 
and 98/iT, respectively, for the two resonances observed 
at 85.3mT (853G) and 90.7mT (907G) |J|. These val- 
ues of A agree with the measured value for the 90.7mT 
resonance |4|, and with the indirectly inferred order of 
magnitude for the 85.3mT resonance Q|. The scatter- 
ing lengths a m and a am for molecule-molecule and atom- 
molecule collisions, respectively, are not known, but cal- 
culations show that the results of our analysis are prac- 
tically insensitive to the variation of their values as long 
as they stay within the order of magnitude of a a . 

In the calculations one should make a clear distinction 
between the two types of experiments conducted at MIT 
- the slow-sweep and the fast-sweep experiments. In 
the first case, the values of the magnetic field at which 
the ramp stops short of resonance are such that the con- 
ditions needed for the excitation mechanism to occur do 
not exist, and only the collisional deactivation applies. 

The graphs shown in Figs. |l| and ^ pertain to the slow- 
sweep MIT experiment for the strong 90.7 mT resonance. 
This resonance has been approached from below with two 
ramp speeds, and from above with one. Figure |l| shows 
the surviving atomic density n, and Fig. |] shows K3, 
both vs. the stopping value of the magnetic field B. 
The difference between Eqs. (|3^), (|4p| ) and the results 
of a direct numerical solution of Eqs. ( |32| ) for all ramp 
speeds is so small that the corresponding plots are indis- 
tinguishable. The calculated plots were obtained using 
homogeneous-density initial conditions, starting from a 
B value of 89.4mT on approach from below and 91.6mT 
from above. The corresponding initial mean densities 
were extracted from the experimental data Q. A best 
fit of the parameter 7, using Eq. (fh)]), to the MIT data 
gives 7 = 0.8 x 10~ 10 cm 3 /s (see Fig. 0). But owing 
to the large scattering of the experimental data, and the 
associated uncertainty in the value of 7, we proceeded us- 
ing the value of 10 _10 cm 3 /s in our calculations, following 
Ref. ||. Given a density of about 10 15 cm -3 this value of 
7 implies a deactivation time of ~ 10~ 5 s. The molecule- 
molecule deactivating collisions are negligible compared 
to the atom-molecule collisions due to the small molecu- 
lar density [see Eq. (|6|)] whenever the fast decay condi- 
tion @ holds. 

The rate of condensate loss due to atom-molecule de- 
activating collisions was also studied in Ref. || . As noted 
in the previous section, the expression for K3 obtained 
there is smaller by a factor 1.5 from that of Eq. ([39]) [see 
discussion after this equation]. Without this factor one 
would not obtain the almost perfect match between the 
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analytical and the numerical results mentioned in the dis- 
cussion of Figs. |] and || above. This omission has been 
corrected in Ref. || and the value they obtain for their 
G s t a b (corresponding to our 27), of 4 x 10 -10 cm 3 /s, was 
estimated by best agreement with the experimental data 
|| . This value is 2.5 times bigger than our estimate. This 
discrepancy may be due to the large scatter in the exper- 
imental data. The estimate of Ref. JB[ is based only on 
the experimental K3 plot (see Fig. 0) which is obtained 
by a differentiation of the experimental density data, and 
therefore shows a scatter of the data of up to an order of 
magnitude (see Fig. |[ as well as the corresponding figure 
in Ref. Q, that use a logarithmic scale). Equation J4C| ) 
(presented in Ref. allows us to estimate the value 
of 7 directly from the experimental plot for the atomic 
density (see Fig. |[), that shows a much smaller (w 20%) 
scatter of the data points. Anyhow, in both cases, the 
error bar should be comparable to the value of 7 itself. 

An inelastic rate coefficient 27 with a magnitude of 
the order of 10 -10 cm 3 /s appears to be reasonable. First, 
this value is two orders of magnitude smaller than the 
upper bound set by the unitarity constraint on the S- 
matrix [ fL9| |. In the limit of small momentum, unitar- 
ity provides 27 < hX/m, where A (the de Broglie wave- 
length) in the current situation is limited by the exper- 
imental trap dimensions. This constraint sets an upper 
bound of 2.5 x 10 _8 cm 3 /s to 27. Second, our estimate of 
10~ 10 cm 3 /s for 7 is consistent with the order of magni- 
tude of recently calculated vibrational deactivation 
rate coefficients due to ultracold collisions of He with H2 
in highly excited vibrational levels. 

The remaining figures (@to [?]) pertain to the fast-sweep 
MIT experiment. Figures^] and ^ present the surviving 
part of the trap population after passing through each 
of the two resonances. Following the experimental condi- 
tions, the value of 10 15 cm -3 is used for the initial density, 
and the magnetic field starting and stopping values are 
shifted from resonance by 3 mT for the strong (90.7 mT) 
resonance and by 2 mT for the weak (85.3 mT) one. 

The analytical results of Eq. ( ]43f), together with the 
direct numerical solutions of EqsT(p2[) for the homoge- 
neous initial distribution, considering only the deactivat- 
ing mechanism, are compared in Fig. |^ with the results 
of the fast-sweep experiment [||4| . Although Fig. ^| does 
not specify the direction of the Zceman shift fiB, one 
should recall (see Sees. || and II C 2 above) that the ex- 
citation mechanism can be ignored when \iB < 0. The 
asymptotic result (^3|) reproduces a characteristic depen- 
dence on the ramp speed, although it is independent of 7. 
The numerical solutions clearly show a dependence on 7, 
as assumptions ([35]) and (Ell) underlying the asymptotic 
result ( p3| ) do not hold in this case. The loss reaches 
a maximum at a value of 7 dependent on the various 
parameters (e.g., 7 « 10~ 10 cm 3 /s for the 90.7 mT reso- 
nance when sno sa 2) , as a result of the conflicting asymp- 
totic and fast-decay conditions ([35]) and The cal- 
culated drop in the condensate loss on increase of the 



molecular dumping rate j m , which may seem paradox- 
ical, can have the following explanation. Reaction (||) 
leads to the loss of three condensate atoms per each reso- 
nant molecule formed, while reaction (||) leads to the loss 
of only two condensate atoms. In the present case, the 
loss rate is limited by resonant molecule formation [re- 
action (0)]. Therefore, the increase of 7 m transfers flow 
from channel (||) to channel (j^), thus reducing the num- 
ber of condensate atoms lost per each resonant molecule 
formed. 

The results of the calculations, in which the effect of 
crossing to excited trap states is incorporated, are plotted 
in Fig. ||. The crossing rate was calculated with Eq. p5|), 
using the values of A, a a , and /1 given at the beginning 
of this Section. The results clearly show that the two 
mechanisms do not augment each other. Adding the two- 
body excitation mechanism to the more efficient three- 
body deactivation mechanism actually reduces the loss. 
A possible explanation of this paradoxical result is similar 
to the one used in explaining Fig. [j|. 

Our predicted losses are somewhat higher than the 
ones obtained in the experiments for the 85.3 mT res- 
onance, but significantly lower for the 90.7 mT reso- 
nance. Actually, the closest we can get to the experi- 
mental results is by considering only the two-body ex- 
citation mechanism for the 85.3 mT resonance, and the 
three-body deactivation mechanism for the 90.7 mT res- 
onance. 

The results of Ref. Q, considering only the two-body 
mechanism, also show a better agreement for the 85.3 mT 
resonance. There are, however, significant differences be- 
tween the theory used there and the one presented here 
and in Ref. jl(J. In Ref. || the parameter 70 (analogous 
to our r cr ) is considered as a constant, independent of 
the released energy E or time t, and prevailing all along 
the sweep, below and above the resonance. The time in- 
tegral in their Eq. (4) should be smaller by a factor of 2 if 
the correct energy dependence of 70 is used. In our ver- 
sion, following Ref. |L(J, r cr is t -dependent, exists only 
for E > 0, and attains the correct Wigner limit when 
E — ► 0. The final result of Ref. || is an expression simi- 
lar to our Eq. ( [43| ) , with the exception that our parameter 
s is larger by a factor 3/2 than the corresponding param- 
eter in Ref. . This difference reflects the fact that Eq. 
d43| ) was derived for the three-body deactivation process, 
while Ref. deals exclusively with a two-body excita- 
tion process. As a matter of fact, their coefficient should 
be further reduced by a factor of 2 associated with the 
energy dependence discussed above. Therefore the mech- 
anism of Ref. H requires a factor s 3 times smaller than 
the one used in Eq. ( ]43| ) for the deactivation mechanism. 

Our results also differ from those of Ref. The 
present calculations take into account the decrease of the 
condensate density during the crossing [see discussion af- 
ter Eq. d43|)] and therefore produce a lower condensate 
loss than the one obtained in Ref. |L0| , in which the loss 
is described by a Landau-Zener formula. In the limit of a 
small loss, Eq. (E3) may resemble a Landau-Zener expres- 



9 



sion in which sn is substituted for the exponent. How- 
ever, a study of the associated free-atom problem shows 
that, even when we retain the initial density n (r, — oo) 
in the exponent, the latter would still be smaller than 
sn by a factor of 3. This factor is directly related to the 
many-body character of the Gross-Pitaevskii equations. 

At the first stage of the atomic condensate loss ([!]) a 
condensate of molecules in the resonant state is formed. 
This molecular condensate is unstable, due to the deac- 
tivation of the molecules by reactions (g) and (||) , as well 
as their decay through the excitation mechanism. Figure 
H presents the calculated time dependence of the molec- 
ular condensate density, for various magnetic field ramp 
speeds, taking into account the various loss mechanisms. 
The oscillations of the molecular condensate density are 
connected to the intercondensate tunneling considered 
in Ref . || . Figure || shows that the excitation loss en- 
hances the damping of the oscillations and the decay of 
the molecular density. Nevertheless, the molecular con- 
densate persists at least a few tenths of a microsecond 
before decaying. This time is long enough to allow con- 
verting the population of the vibrationally-excited state 
to to the ground molecular state by methods of coherent 
control |2(j,^l|. Various techniques exist today for trans- 
fering populations coherently to a preselected state p|| . 
The choice of technique would be dictated by properties 
of the process and of the target state (such as selection 
rules and stability). 

The calculated peak values of the resonant molecular 
state density are presented in Figs. || and R for various 
rates of the deactivating collisions. Figure o refers only 
to the deactivation mechanism involving atom-molecule 
and molecule-molecule collisions, which take place for 
fiB < 0, whereas Fig. takes into consideration the com- 
bined effect of the two mechanisms (deactivation and ex- 
citation) which may take place for ijlB < 0. These figures 
show that from about 10% to 90% of the atomic density 
can be converted temporarily to a molecular condensate 
in the resonant state, in spite of losses due to the excita- 
tion mechanism and the molecule-molecule deactivation 
collisions. The calculated molecular condensate density 
is higher for the 90.7 mT resonance, or when the mag- 
netic field ramp speed is not too large. 

V. CONCLUSIONS 

This paper discusses the two types of loss experiments 
HQ conducted at MIT on sodium BEC, using a time- 
varying magnetic field in the proximity of a Feshbach 
resonance. The various processes discussed here involve 
a temporary formation of a molecular condensate. In the 
slow-sweep experiment the dominant loss mechanism is 
three-body deactivation by atom-molecule inelastic colli- 
sions. A best fit of an analytical expression for the atomic 
density obtained here [Eq. (^p|)1 to the MIT data yields a 
rate coefficient for deactivating atom-molecule collisions 



of 2j — 1.6 x 10~ 10 cm 3 /s. For the fast-sweep experi- 
ment an analytical expression, generalizing the Landau- 
Zener formula to a case of coupled nonlinear equations, 
is obtained. A different two-body mechanism, involving 
an excitation of the condensate by curve crossing, has 
been proposed previously JlTj] . The combined effect of the 
two mechanisms is studied here, including also molecule- 
molecule deactivating collisions. The analysis shows that 
both processes should be taken into account, that they 
do not contribute additively to the loss, and that the 
outcome of the competition between them varies from 
one Feshbach resonance to another. Our numerical re- 
sults show that, under favorable conditions, a substantial 
fraction of the trap population is converted to an unsta- 
ble molecular condensate. This condensate persists long 
enough to allow its coherent transfer to a more stable 
state. 
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.4 90.7 

B(mT) 

FIG. 1. The surviving mean density vs. the stopping value 
of the magnetic field, calculated with the optimal value of the 
deactivation parameter 7 = 0.8 x 10~ n cm 3 /s. The reso- 
nance was approached from below with two ramp speeds, 13 
mT/s and 31 mT/s, or from above, with the ramp speed -6 
mT/s. These are compared with the experimental results Q] 
(squares, triangles and circles) for several values of the ramp 
speed dB/dt (in mT/s). 

1 0= 




90.7 
B(mT) 

FIG. 2. The 3-body rate coefficient (K3) vs. the stop- 
ping value of the magnetic field, calculated with the optimal 
and other values of the deactivation parameter 7 (in units of 
cm 3 /s), on approaching the resonance from below or above. 
The calculated value of Kz is independent of the ramp speed 
value. Other notations as in Fig. 1. 
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FIG. 3. Ratio of surviving trap population N to the initial 
one iV for the 85.3mT (853G) and 90.7mT (907G) resonances 
[parts (a) and (b), respectively] in the homogeneous-density 
approximation vs. sno (where the parameter s is defined by 
Eq. ([I3]) and no is the initial density) . The curves show results 
of calculations carried out for different magnitudes of the pa- 
rameters 7 and 7 m (in units of cm 3 /s), without accounting for 
the excitation mechanism. The asymptotic analytical result 
Eq. (|5||) is given by the dotted line. The results of the MIT 
fast-sweep experiment (Ref. 14]) are shown for comparison. 
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FIG. 4. Same as Fig. |3| but accounting for the excitation 
mechanism. The results of calculations without account- 
ing for the excitation mechanism are given for comparison 
by the dotted line. In the part (b) the plot calculated for 
7 = 10 _10 cm 3 /s and different magnitudes of j m are practi- 
cally indistinguishable. 
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FIG. 5. Time-dependence of the molecular condensate den- 
sity in cm -3 during fast passage through the 85.3mT (853G) 
and 90.7mT (907G) resonances [parts (a) and (b), respec- 
tively], calculated using the value of 7 = 10~ 10 cm 3 /s. The 
dashed and solid lines are calculated accounting for the com- 
bined effect of the two mechanisms (deactivation and excita- 
tion) for sno = 1 and 5, respectively, and 7 m = 10 _9 cm 3 /s. 
The dotted and dash-dotted lines refer only to the deactiva- 
tion mechanisms for 7 m = and 10 _9 cm 3 /s, respectively, and 
sno = 1. The oscillations of the molecular density are faster 
for the lower value of sno. 
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FIG. 6. Peak values of the molecular condensate density 
in cm -3 attained during fast passage through the 85.3mT 
(853G) and 90.7mT (907G) resonances [parts (a) and (b), re- 
spectively] vs. sno for various magnitudes of the parameters7 
and 7 m (see Caption of Fig. ^|), calculated without accounting 
for the excitation mechanism. 
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FIG. 7. Same as Fig. but accounting for the excitation 
mechanism. The results of calculations without account- 
ing for the excitation mechanism are given for comparison 
by the dotted line. In the part (b) the plot calculated for 
7 = 10~ 10 cm 3 /s with 7 m = and 10~ 10 cm 3 /s are practi- 
cally indistinguishable. 



TABLE I. Asymptotic conditions as consistency tests on solutions of Eq. ( |37| ) 





case 


h 


n 



a n -^const h — > n -^const 

b n/t -» h ~ -6|g| 2 7n 3 / (/iBt) 7i — > const 

c n/t — > oo n ~ — 6|g| 2 n/ (^7) n ~ exp [— 6|<?| 2 t/ (?i 2 7)] 

d n/t -^const h ~ t n ~ t 2 
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